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WHOLE-BEAM  THERMAL  BLOOMING  EXPERIMENTS  WITH  LONG  LASER  PULSES 

Liu  Changling,  Hu  Xijing,  Yuan  Yonghua . 

Gui  Yuanzhen,  Zhang  Kexing,  Li  Qiming, 

Wang  Chunyan,  Liu  Xufa,  and  Yang  Libing 

Southwest  Institute  of  Fluid  Physics 
P.O.  Box  523,  Chengdu  610003 

ABSTRACT  Whole- beam  thermal  blooming  experiments  have  been  done  using  a  focused  long 
pulse  laser  (t^“  1 .8ms  .  £  =  200J»  angle  spread  (?  =  2mrad  ) .  A  blooming  cell  is  full  of  ammonia  at 
0^6  atmosphere  pressure.  A  CCD  system  is  used  to  measure  the  diameter  of  the  distorted  laser 
beam  waist  and  the  intensity  profile  of  the  spot  at  waist.  The  absorption  cofficients  are  measured 
by  a  digital  calorimeter.  The  instantaneous  irradiation  intensity  at  spot  centre  is  measured  by  a 
photodiode  digital  collector  -  Roland  mapper-Computcr  system.  The  instantaneous  intensity  pro¬ 
file  is  measured  by  a  high  speed  camera. 

KKV  WORDS  absorption,  blooming,  distortion.  laser  atmosphere  propagation. 


I.  Introduction 

The  reason  for  paying  attention  to  the  transmission  of  high- 
powered  lasers  in  the  atmosphere  is  mainly  because  these  can 
become  a  new-generation  strategic  weapon.  One  of  the  main 
approaches  in  attacking  targets  in  space  is  to  induce  structural 
failure  of  the  target  with  the  thermal  and  mechanical  effects  of 
the  laser.  High-powered  laser  systems  are  required  to  produce 
such  destruction.  Therefore  the  authors  are  studying  the 
nonlinear  effects  of  high-powered  lasers  in  atmospheric 
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transmission;  there  is  practical  significance  such  as  absorption, 
scattering,  atmospheric  turbulence,  thermal  blooming,  gas 
puncturing,  excited  Raman  scattering,  and  optical  compensation. 

The  effect  of  thermal  blooming  is  due  to  absorption  of  laser 
energy  in  the  laser  beam  by  the  atmosphere  so  that  the  initial 
medium  in  thermal  equilibrium  is  heated  locally  to  cause  an 
increase  in  gas  density  and  raise  the  speed  of  sound,  thus 
inducing  changes  in  refractivity  and  causing  a  negative-lens 
effect  to  scatter  the  laser  beam.  Thus,  the  light  beam 
dimensions  are  increased  at  the  focal  plane  to  cause  light-spot 
anomalies.  For  systematic  research,  frequently  some  absorptive 
media  are  used  to  simulate  atmospheric  absorption  for  thermal 
blooming  research  in  laboratories.  There  are  different  kinds  of 
thermal  blooming  with  different  laser  acting  periods.  When 
classified  by  time  properties,  there  are  continuous-wave  lasers, 
pulsed  lasers,  and  repetition-frequency  pulsed  lasers.  There  are 
two  thermal  blooming  situations  for  single-pulse  lasers.  Which 
situation  is  involved  is  determined  by  the  time  (of  flow  dynamic 
properties)  Tj.=D/Cj  .  D  is  the  light  beam  diameter;  Cj  is  the 
speed  of  sound  in  the  medium;  and  Tjj  is  the  time  required  for 
sound  wave  to  pass  through  the  light  beam  cross-section  in  the 
absorbing  medium.  When  the  laser  pulse  time  tp<T[j,  the  density 
disturbances  in  the  laser  beam  increases  as  the  third  power  of 
time;  this  is  called  the  tertiary  thermal  blooming  effect.  When 
tp>t[j,  the  density  disturbances  increase  in  proportion  to  time. 

In  this  experiment,  what  decided  by  the  light-beam  diameter. 


2 


Tj,<<tp  is  the  first-power  of  the  thermal  blooming  effect. 

Detailed  experimental  and  theoretical  studies  were  conducted 
between  the  thermal  blooming  effect  and  time  in  the  transmission 
of  high  power  lasers  that  were  carried  out  abroad.  In  [1,2],  a 
CO2  laser  operating  at  10 . 6micrometers  was  used.  In  references 
[3,4] ,  measurements  of  the  time  process  of  a  1 . 06micrometer 
neodymium  glass  laser  were  carried  out,  thus  yielding  the 
relatively  consistent  results  compared  with  the  theoretical 
estimates.  In  1973,  at  Lincoln  Laboratories  [4]  in  the  United 
States,  the  1 . 06micrometer  pulsed  laser  (Ep=20  to  200J,  |3=50)  , 

NHj  of  O.SMPa  was  used  as  the  absorbing  medium  to  simulate 
absorption  in  the  atmosphere.  Experiments  on  the  time  process  of 
various  laser  pulses  were  conducted.  In  the  authors' 
experiments,  thermal  blooming  experiments  and  research  were 
conducted  by  using  1.8ms  of  pulse  time,  200 J  of  output  energy,  in 
the  f ree-oscillation  neodymium  glass  laser  of  light-beam 
divergence  angle  of  2mrad  ((3=20),  in  free  oscillations  with  lens 
focusing  while  passing  a  thermal  blooming  cell  filled  with 
ammonia  at  0  to  O.SMPa.  A  CCD  system  was  used  to  measure  the 
variational  curve  between  the  time-integrated  light  meter  and  the 
ammonia  pressure  of  the  thermal  blooming  cell  of  the  laser  beam 
waist.  With  an  energy  meter  that  had  a  digital  display  and  peak 
value  retention,  the  absorption  coefficient  was  measured  after 
the  laser  beam  had  passed  through  the  ammonia  cell.  A  complex  of 
silicon  optoelectronic  diode/digital  collection  and  analysis 
instrument/Roland  automatic  plotter/computer  processing  system 
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was  used  to  determine  the  time-dependent  variation  curve  and 
light-intensity-gas  pressure  in  amm.onia  cell-time  curve,  for  the 
central  light  intensity  of  the  laser. 

By  using  the  KDP  frequency  doubler/high— frequency  frame 
dividing  camera  system  to  record  the  light  intensity 
distribution— time  image  after  the  laser  beam  had  passed  through 
the  ammonia  cell,  thus  we  obtained  frame-by-frame  photographs, 
and  observation  of  the  time-dependent  anomaly  of  the  light-beam 
field  intensity.  Thus,  beam  broadening  and  central  hollowing 
phenomenon  (the  intensity  in  the  light  beam  center  decreased,  and 
the  intensity  at  the  beam  periphery  was  greater  than  the  center 
intensity) ,  after  the  light  spot  had  passed  through  the  thermal 
blooming  cell.  The  experimental  results  are  basically  in 
agreement  with  theory. 


II.  Brief  Description  of  Analytical  Theory  [5] 

Based  on  the  conservation  of  energy,  in  the  approximate 
conditions  of  the  side  axis,  the  light  intensity  I  satisfies  the 
following  equation: 


V,nd:V,I  + 


V  nd:  V  I  ^ 


I 

II 


V  ;  II  dz=  -  il 


(  I  ) 


In  the  equation,  n  indicates  the  refractivity  of  the  medium;  r 
indicates  the  horizontal  direction;  and  p  indicates  the  total 
absorption  coefficient  of  the  medium. 

Under  the  conditions  of  neglecting  the  longitudinal-gradient 
of  the  refractive  index,  neglecting  the  diffraction  effects, 
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(that  is,  the  application  of  geometric  optics),  Eq.  (1)  can  be 
written  in  integrated  form 


I[r,z) 

/(r.O) 


=  exp{ 


iz  - 


I'x  -V,/  /) 


t  n 


iiz'dz' ; 


(  2  ) 


The  variations  of  the  refractive  index  n  and  density  p  of 
the  gas  with  time  satisfies  the  Dale-Gladesrone  law: 


«  -  1  =  X’  /' 


t  3  ) 


In  the  equation,  k  is  the  proportional  constant. 

In  the  transmission  problem  of  single-pulse  lasers,  it  is 
required  to  consider  the  variation  of  density  with  time.  Based 
on  the  equation  sets  in  fluid  mechanics,  the  density  p  satisfies 

(4, 

In  the  equation,  c^  is  the  speed  of  sound  in  the  absorbing 
medium;  f) .  c,  ,  c^  is  the  ratio  of  specific  heat;  and 

p  is  the  gas  pressure  of  the  absorbing  medium. 

The  equation  for  the  ideal  state  of  the  gas  is 

p  =  pRT  .  R  =  i  J\ )  (  5  ) 


The  problems  being  studied  by  the  authors  include  the  fact 
that  the  first  term  tp>tj|  in  Eq.  (4)  can  be  neglected.  By 
integrating  Eq.  (4)  with  respect  to  time,  we  have 


=  - 


7-1  1 

l(U= - - - 7—  %  - 


71  a 


R(  t )  (It 


If  a  laser  beam  is  a  rectangular  pulse,  then  Eq.  (6)  is 
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With  respect  to  a  gaussian  light  beam,  from  Eqs .  (2),  (3), 

and  (6),  we  can  obtain  the  analytical  solution  for  light- 
intensity  distribution  after  passing  through  the  absorbing 
medium : 


/(r.r.r) 
/(r.O  ,  O 


=  e““‘exp 


•^0 
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(5)  , 
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In  the  equation 


5(ctz)=  —  [1-  —  (l-e-“)] 
az  az 


The  anomaly  parameter  of  defining  the  alignment  beam 


2n^a  P{t)dtz^ 
- - -  f(az) 


Eq.  (8)  is  transformed  to  become 
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In  Eq .  (10)  ,  nj  is  dn/dT;  this  is  a  negative  value  for  nip  of  most 

gases,  and  p  is  the  absorption  coefficient,  P{t)  is  the  laser 
power  of  the  laser  device;  z  is  the  length  of  the  absorbing 
medium;  ng  is  the  refractive  index  of  the  gas  without 
disturbances;  pg  is  the  initial  density  of  the  absorbing  medium; 
Cp  is  the  isobaric  specific  heat.  p  is  the  beam  radius  when  the 
light  intensity  is  reduced  to  1/e. 

In  Eq .  (11),  l(r,0,t)  is  the  intensity  distribution  of  the 

laser  beam  without  thermal  blooming.  For  a  gaussian  light  beam, 
we  have 


(11  ) 
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In  the  equation,  Po(t)  is  the  initial  power  of  the  laser  beam. 

Eq.  (11)  is  the  theoretical  result  of  the  medium  with  respect  to 
the  straightening  situation.  This  should  be  revised  for  a 
focusing  beam.  Assume  that  the  focusing  beam  satisfies  the 
geometric-optics  approximation,  in  other  words. 


(13) 


Refer  to  Fig.  1  for  the  light-path  diagram.  In  the  equation,  ap 
is  the  light  beam  initial  radius,  a^  is  the  radius  of  the  focal 
spot;  Zf  is  the  focal  distance;  and  a^  is  the  beam  radius  at 
distance  z. 


For  long  pulses  (tp>Tjj) 


(14) 
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,  ,  x^-l-2lnx  a 

^{.r)=  - ^ -  .  ;,=  iL 

2(1  -  — 

X 

When  x=l,  q=l .  When  x>>l,  q{x)  is  proportional  to  X2/2. 

In  Eq .  (13) ,  a^  is  the  beam  radius  at  z=0 .  aj  is  the  beam  radius 

at  z=L  (exit  of  blooming  cell). 

The  temperature  rise  generated  by  the  medium  when  absorbing 
laser  energy  is 


(15) 


AT- 


a£ 


( 16) 


In  the  equation,  the  unit  of  «  is  m'^-  E(J)  is  the  pulse  energy; 

2 

and  S  (m  )  is  the  cross-sectional  area  of  the  light  beam.  The 

rate  and  time  of  thermal  expansion  [5]  are 

K  D 

y  ^  - 

P  C.  CJ  V  , 

K  (W  m  K^)  is  the  coefficient  of  thermal  expansion.  The 
convection  rate  and  the  property  time  induced  by  the  floating 


force  are  given  by 


C=( 


2a  Pg 

pc,T 


r 


(IS) 


In  the  equation,  g  is  acceleration  due  to  gravity  and  D  is  the 
beam,  diameter. 

The  characteristic  time  in  fluid  mechanics  is  given  by 


From  Eq .  (13)  we  obtain: 


r 


n 


D_ 


a,  =  0.867cm  ,  =  502cm. 


(19) 


From  Eq.  (15) ,  the  correction  factor  (of  the  focusing  beam) 
q (x) =2 . 508 . 
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III.  Whole  Beam  Thermal  Blooming  Experimental  System  of  a  Long- 
Pulse  High-Powered  Laser 

From  Fig.  2,  the  overall  layout  of  the  experimental  system, 
we  can  see  that  a  free-oscillating  neodymium  glass  laser  device 


Fig.  2  Whole- beam  blooming  experimental  layout 


with  200J  output  energy  and  1.8ms  pulse  time  is  used  in  order  to 
carry  out  the  experimental  research  on  whole-beam  thermal 
blooming  and  in  obtaining  the  required  parameters  of  the  thermal- 
blooming  effect.  A  blooming  cell  with  ammonia  at  the  highest 
pressure  0 . 6MPa  was  used  with  the  optimal  beam  diameter  at  OD20mm 
and  optical  beam  divergence  angle  of  2mrad  for  the  output  of  a 
parallel  plane  cavity.  The  blooming  cell  consists  of  an 
ODD . Imxlm  stainless  steel  tube  flange  system,  valve,  and  pressure 
gauge  for  measuring  the  ammonia  gas  pressure  in  the  cell.  Both 
ends  of  the  tube  are  sealed  with  optical  windows  of  high  optical 
quality.  The  ammonia  gas  in  the  cell  is  used  to  simulate  the 
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laser  of  the  absorbing  inediuin  of  the  atmosphere.  This 
arrangement  is  used  to  focus  the  high-powered  laser  beam  after 
focusing  and  passing  through  the  2m  focal  length  of  the  blooming 
cell  to  focus  onto  the  convergent  lens  Lj  and  other  transmission 
and  reflective  optical  components. 

With  respect  to  the  test  system  after  the  focusing  light 
beam  passes  through  the  light  cell,  the  test  system  is  used  for 
the  variation  of  ammonia  pressure  in  the  blooming  cell  with  the 
time-integrated  diagram  of  the  laser  beam  waist.  The  absorption 
coefficient  test  system  for  variation  in  ammonia  pressure  after 
the  focused  laser  beam  has  passed  through  the  blooming  cell.  The 
following  four  components  constitute  a  controlled  system 
synchronized  with  the  experiment:  the  central  intensity  of  the 
laser,  with  dependence  on  the  laser  pulse,  measurement  processing 
system  varying  with  ammonia  gas  pressure,  distribution  of  field 
intensity  of  the  laser  beam  cross-section  after  passing  through 
the  blooming  cell,  and  the  test  system  for  development  variation 
of  anomalies  with  time.  After  the  laser  is  focused  through  the 
convergent  lens  and  passes  through  the  blooming  cell,  when  there 
is  no  ammonia  in  the  cell,  the  output  of  the  light  spots  with 
relatively  concentrated  field  intensity  is  shown  with 
photographic  paper  recording  the  laser  beam  waist  (Fig.  3).  When 
the  cell  is  filled  with  0.1  to  O.GMPa  ammonia  gas  from  low  to 
high,  these  light  spots  gradually  diverge,  decay,  and  are 
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anomalous  with  respect  to  pressure  intensity  (Figs.  3b,  c,  d,  e, 
f ,  g)  . 


1  lu-  Spot.N  vM  l.iNcr  ujui  on  phoh»iiraphjc’  prinJing  paper  lor  t)  O.oMPa  anmuMita 
pressure  in  the  blooming  cell 


IV.  Experimental  Tests  and  Analysis  of  Thermal  Blooming  Effects 
4.1.  Time— integrated  Image  of  Focused  Laser  Beam  Waist 

As  shown  in  Fig.  1,  after  reflection  with  convergence 
through  lens  Lj ,  the  200J  pulsed  laser  beam  passes  through  an 
optical  wedge  M,  with  single-plane  reflectivity  of  4%,  reflective 
mirrors  Mj  {R=95%)  ,  and  Mj  (R=95%)  ,  the  beam  waist  of  the  light 
beam,  falls  onto  an  unfinished  glass  screen  in  front  of  the 
optical  lens  of  the  CCD  camera.  This  laser  beam  does  not  pass 
through  the  blooming  cell.  The  portion  of  the  laser  beam  passing 
through  the  optical  wedge  Mj  and  after  passing  through  a  blooming 
cell  filled  with  ammonia  gas  at  a  different  pressure  and 
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reflected  through  optical  wedge  {R=4%) ,  the  beam  waist  of  this 
laser  beam  also  is  incident  on  the  unfinished  glass.  The  light 
intensity  distribution  diagram  of  both  laser  beam  waists  are 
integrated  on  the  CCD  during  the  entire  time  period  of  1.8ms,  and 
is  collected  and  processed  with  a  computer  connected  to  the  CCD 
camera.  Thus,  on  the  computer  screen  the  color  two-  and  three- 
dimensional  field-intensity  distribution  diagrams  of  time 
integration  of  two  laser  beam  waists,  one  passing  through,  and 
the  other  not  passing  through  the  blooming  cell;  the  diameters  of 
the  two  laser  beam  waists  are  shown.  In  Fig.  4a  and  b  are  shown 


( )  (  h  ) 

!  linic"  nUcrjjrjlcd  -D  jnlcn.Mlv  .11  hcani  v^aist  mca.'iurcd  by  CCl)  sy.sicni 


two  dimensional  color  light  intensity  distribution  diagrams  and 
the  derived  data  of  beam  waist  diameters  for  the  laser  beam  waist 
of  the  blooming  cell  filled  with  0.1  to  0 . 5MPa  ammonia.  When  the 
blooming  cell  is  filled  with  0  to  O.SMPa  ammonia,  respectively, 
then  the  field  intensity  distribution  diagrams  of  seven  beam 
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waists  and  diameter  data  for  the  beam  waist  are  obtained  as  shown 
in  Fig.  5  and  Table  1. 


Table  I  Diametcrd  at  beam  waist 
«  1 


P  O.IMPa 

U 

I 

3 

4 

5 

6 

d  mm 

4 

4.2S 

4,5! 

6.S8 

8.83 

9.83 

13 

Kig.  Curve  ol'diameler  at  beam  waist  versus 
ammonia  pressure 


4.2.  Laser  Beam  Absorption  Coefficient  When  Beam  Passes  Through 
Bloom.ing  Cell 

From  Fig.  2  the  laser  beam  not  passing  through  the  blooming 
cell  is  reflected  by  Mj  and  ,  and  then  passes  through 
(T=0.05)  and  Mj  (T=0.92)  to  arrive  at  the  energy  meter  Ej .  The 
other  laser  beam  passes  through  the  front  window  opening  (T=0.92) 
of  the  blooming  cell  and  {T=0.92) ,  then  passes  through  the 
blooming  cell  before  passing  through  the  rear  window  opening  of 
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the  blooining  cell,  wedge  lengths  Mj  and  M,  {T=0.92),  and  then 
reflected  (R=0.04)  from  the  front  surface  of  the  wedge  length 
to  arrive  at  the  energy  meters  Ej  and  Ej .  Records  of  the  energy 
recording  instruments  are  stored  as  digital  peak  values  with  E2. 
Laser  energy  Ej  not  passing  through  the  blooming  cell,  and  Ej 
passing  through  the  blooming  cell  are  measured,  respectively. 
Since  both  light  beams  pass  through  different  optical  systems,  it 
is  required  to  divide  with  the  transmission  coefficients  the 
measured  energy  Ej  and  E2  of  both  energy  meters  to  obtain  the 
actual  energy.  Since  there  are  system  errors  in  the  two  energy 
meter  systems,  it  is  required  to  calibrate  {in  order  to  remove 
system  errors)  of  the  two  energy  test  systems  before  the  blooming 
cell  is  filled  with  ammonia  and  in  the  condition  of  laser  device 
output  of  a  certain  amount  of  energy.  The  blooming  cells  are 
filled,  respectively,  with  ammonia  gas  at  0  to  O.SMPa.  Under  the 
conditions  of  200J  energy  outputted  from  the  laser  device,  the 
absorption  coefficients  (Table  2)  are  obtained  for  passage 

Table  2  Absofbtion  coefficients  in  blooming  cell 


WMPa 

0 

O.i 

0.2 

0.3 

0.4 

0.5 

0.6 

2/  m  ■' 

0.0054 

0.076 

0.21 

0.54 

0.68 

0.88 

1.06 

through  the  Im  blooming  cell.  Thus,  the  curve  of  the  absorption 
coefficients  was  obtained  for  a  laser  by  the  ammonia  blooming 
cell  of  0  to  O.SMPa  (Fig.  6). 
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4.3.  Time  Variation  Curve  of  Central  Light  Intensity  of  Laser 
Beam 

The  time  variations  of  laser  pulses  versus  central  light 
intensity  (Fig.  7)  of  the  laser  after  passing  through  the  0  to 
0.6MPa  ammonia  blooming  cell  are  recorded,  respectively,  with  a 
silicon  light-emitting  diode/digital  collection  and  analysis 


l  iy.b  Curve  oC  ;ib5orption  collicicm  wr-ius 
ammonia  pressure 


instrument/automatic  plotter  (Fig.  2) .  From  the  wave  curves  of 
these  six  lasers  we  can  see  that  absorption  is  higher  for  a  laser 
beam  passing  through  a  cell  with  higher  ammonia  gas  pressure. 
After  the  laser  beam  passes  through  the  blooming  cell  filled  with 
0.3MPa  gas,  the  central  light  intensity  of  the  laser  of  1.8ms 
laser  pulse  time  has  its  peak  intensity  reduced  to  approximately 
1.3  at  the  time  of  400microseconds .  After  1.2ms,  the  intensity 
drops  to  almost  zero.  When  passing  through  the  blooming  cell 
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Fig.  7  Curve  of  intensity  at  the  centre  versus  lime,  the  values  across  curves  indie. tic 
atmosphere  pressure  inside  blooming  cell 

filled  with  ammonia  at  a  gas  pressure  of  0.6MPa,  the  central 
light  intensity  of  the  laser  beam  drops  to  one-third  of  the  peak 
value  intensity  only  after  200microseconds .  After  1ms,  it  drops 
to  almost  zero.  The  light  intensity  curve  I  and  time  t  of  the 
central  time  intensity  curves  of  the  six  light  beams  are 
processed  with  a  computer  system  to  compare  with  the  light 
intensity  Iq  of  the  blooming  cell  with  ammonia  gas  dropping  to  0 
as  the  light  beam  passes  through,  the  time  versus  variation  curve 
of  the  laser  pulse  with  the  central  light  intensity  of  the  laser 
can  be  obtained  (Fig.  8). 

4.4.  Anomalous  Images  of  Light  Intensity  Distribution  of  Laser 
Beam 

When  a  laser  beam  passes  through  a  blooming  cell  filled  with 
ammonia  gas,  as  the  effect  of  thermal  blooming  induces  the 
absorbing  medium  in  the  optical  path  to  generate  a  negative— lens 
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10 


ij  ms 

Fig.  8  Curve  of  //  /,  at  the  centre  versus  time 


effect,  the  light  intensity  of  the  laser  beam  continuously 
increases  with  development  anomaly  of  the  laser  pulse  type.  To 
obtain  the  real-tim.e  im.age  of  the  anom.aly,  a  high-speed  fram.e-by 
frame  camera  of  25,000  frames  per  second  is  used  to  record  the 
im.ages .  Since  the  laser  wavelength  outputted  from  the  neodymium, 
glass  laser  device  is  1 .  OSmicrom.eters ,  light  exposure  is  not 
obtained  for  the  conventional  high-speed  photographic  negative. 
Therefore,  the  authors  used  the  KDP  crystal  frequency  doubler. 

By  using  a  neodym.ium.  glass  laser  device  with  pulse  tim.e 
1 .  Sm.illiseconds ,  when  the  200J  laser  energy  arrives  at  the  KDP 
crystal  frequency  doubler  by  passing  through  the  optical  path  as 
shown  in  Fig.  1,  the  light  intensity  distribution  outputted  with 
doubled  frequency  is  m.easured  with  a  CCD  light  intensity  analysi 
system.  As  indicated  in  the  results,  the  light  intensity 
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distribution  is  proportional  to  the  light  intensity  distribution 
prior  to  frequency  doubling.  Fig.  9  shows  the  high-speed  frame— 


Fig.  9  The  instanianous  irradialion  insien.siiy  of  llie  spot  at  dilTcrcnt  time 
(the  duration  hciweeti  two  spots  in  succession  is  40;is ) 


by-frame  pictures  with  field  intensity  distribution  of  the  laser 
beam  after  the  laser  beam  passes  through  the  thermal  blooming 
cell  filled  with  ammonia  gas  at  a  pressure  of  0.6MPa.  From  these 
frame-by-frame  pictures,  there  is  only  40microseconds  from  the 
second  frame  to  the  third  frame,  great  anomaly  is  generated  in 
field  intensity  distribution,  thus  a  large  empty  area  appeared  in 
the  central  light  intensity  of  the  light  beam.  Hereafter,  the 
light  beam  rapidly  diverges,  and  the  empty  hole  quickly  expands, 
leading  to  lower  and  lower  local  light  intensities  in  the  two- 
dimensional  distribution  of  the  light  beam. 

V.  Discussion  of  Experimental  Results 

Tables  3  and  4  show  the  parameters  of  ammonia  gas  and  laser 
beams  used  in  the  experiment.  In  the  tables,  K  is  the  thermal 
conductivity  coefficient,  and  D  is  the  initial  diameter  of  the 
laser  beam.  Based  on  Eq.  (14)  and  the  absorption  coefficient 
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Tabic  3  NH,  paf^mcier  in  experiment 


( kg  •  m‘  ) 

'I., 

kg  '•  K  ■) 

1  ,  (nr  s') 

■  A', MW  •  m  ’  •  K'') 

0.759 

1.283  X  10  ‘ 

1,000376 

2.241 

415 

0.21 

Table  4  Parameter  laser  set 


^  J 

ms 

/ '  fim 

01  mm 

0‘  mrud 

n 

200 

1.8 

1.06 

20 

2 

20 

a  measured  in  the  experiment,  the  anomaly  parameter  and  the 
temperature  rise  of  ammonia  at  different  pressures  are  listed  in 
Table  5.  From  Table  5,  when  the  ammonia  gas  pressure  is  higher 
than  O.lMPa,  Ng>10,  and  A>10%.  Therefore,  the  analytical  theory 
is  not  applicable  here. 


Tabic  5  Distorlion  parameter  and  risen  temperature  AT  versus  ammonia  pressure 


/>;().  I  MPa 

0 

1 

2 

3 

4 

5 

6 

0.93 

12.7 

34.2 

78.4 

94.6 

115.5 

132.2 

T!K 

17.8 

225 

421 

802 

807 

871 

899 

With  the  estimate  from  analytical  theory  and  the 
experimental  results  of  the  central  light  intensity  of  laser 
beams,  due  to  limitations  of  the  analytical  theory,  the  authors 
can  discuss  only  one  situation,  that  the  ammonia  gas  pressure  is 
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0.  Fig.  10  indicates 


the  time  spectrum  of  laser  intensity  at  the 


Fig.  10  Irradialion  intensity  at  the  centre  of 
Its  waist  versus  time  after  light  beam  goes 
through  a  cell  with  0  pressure 
1 .  laser  source  ;  2.  experimental  result  ; 

3.  analytical  theory  result 

light  beam  center .  This  indicates  that  the  theoretical  curve  and 
the  experimental  curve  are  in  agreement  with  each  other.  For 
comparison,  similar  curves  (Fig.  11)  are  also  given  for  the 
blooming  cell  at  a  gas  pressure  of  O.lMPa.  Here,  the  anomaly 
parameter  has  attained  a  value  of  about  10.  From  Fig.  11,  the 
difference  between  theoretical  and  experimental  values  are 
relatively  high. 

With  respect  to  the  theoretical  and  experimental  values  of 
time  integration  of  the  laser  beam  waist  diameters,  based  on  the 
analytical  theory  in  the  second  section,  time  integration  is 
conducted  on  the  point  r=0  in  Eq.  (8) ,  the  time  integration 
values  on  light  beam  central  intensity  are  calculated, 
respectively,  in  order  to  compare  with  the  central  intensities  of 
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time  integration  measured  with  the  CCD. 


Kig.  1  I  Irradiation  intensity  at  the  tx'nirc  of 
its  waist  versus  time  after  light  beam  goes 
through  a  ceil  with  one  pressure 
I.  laser  source;  2.  experimental  result  ; 

3.  analytical  theory  result 


The  value  of  time  integration  for  light  beam  central 
intensity  is  Iq.  After  passing  through  the  absorption  cell  with 
ammonia  gas  pressure  at  OMPa,  the  time  integration  value  of  the 
light  beam  central  intensity  is  Ij,  then  the  theoretical  value 
Ij/Iq  is  0.44  and  the  corresponding  experimental  value  is  0.63. 

We  can  see  that  the  theoretical  and  experimental  values  are  in 
relative  agreement  with  each  other. 

VI .  Conclusions 

(1)  The  whole-beam  thermal  blooming  experiment  with  long 
pulses  was  successful.  The  absorption  coefficient  of  ammonia  gas 
for  a  light  beam  of  1 . 06micrometer  wavelength  was  measured.  The 
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beam  broadening  and  anomaly  phenomena  of  the  laser  beam  due  to 
the  effect  of  thermal  blooming  was  observed. 

(2)  The  analytical  theory  of  whole-beam  thermal  blooming  can 
be  adaptable  only  to  weak  anomalies.  When  the  anomaly  parameter 
is  approximately  over  10,  the  numerical  simulation  calculation  of 
nonlinear  optics  should  be  applied  in  order  to  precisely  describe 
the  thermal  blooming  phenomenon.  Based  on  the  authors' 
experimental  data,  the  numerical  simulation  calculations  were 
conducted  by  Li  Youkuan  of  the  Beijing  Institute  of  Applied 
Physics  and  Computational  Mathematics,  obtaining  six  curves 
(Fig.  12)  showing  the  time  variation  of  central  light  intensities 


l  ig.  12  Curve  orinlcnsity  al  the  centre  versus  time 

of  lasers  passing  through  a  blooming  cell  containing  ammonia  gas 
at  a  pressure  of  0.1  to  0.6MPa.  The  abscissa  in  the  figure  shows 
19  equal  divisions  when  the  maximum  pulse  time  value  is  1ms.  The 
vertical  axis  is  the  dimensionless  value  of  I/Iq  of  the  laser 
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central  light  intensities.  The  unified  value  of  5.636  is  1.  The 
time-dependent  anomaly  of  three-dimensional  field  intensity 
distribution  of  laser  beams  when  the  ammonia  gas  pressure  is 
0.3MPa  was  also  obtained  (in  Fig.  13,  the  time  gap  between  each 
of  four  frames  of  figures  is  0.3ms).  From  the  comparison  of  the 


I  hi.l.'  iiMcmmU  <it  laser  beam  ver^u^  lime 


numerical  simulation  calculation  results  and  the  thermal  blooming 
experimental  results,  there  is  a  relatively  good  consistency. 

(3)  From  the  results  of  testing  the  absorption  coefficients 
in  Table  2,  when  the  laser  energy  is  about  200J,  the  absorption 
coefficient  is  0.076M^  when  the  gas  pressure  in  the  blooming 
cell  is  O.lMPa.  As  indicated  by  repeated  experiments,  the  data 
obtained  by  the  authors  are  reliable. 


23 


The  first  draft  was  received  on  February  9,  1993;  the  final 
revised  draft  was  received  for  publication  on  May  6,  1993. 


REFERENCES 

[1]  Kenemuth  JR,  Hogge  CB  and  Avizonis  PV.  Appl  Phys  Leii>  1970.  17  (5  ):  220. 

[2]  Buser  RG  and  Rohde  RS,  Appl  Opi^  1975.  14  (  i  ) :  50. 

[3]  Longaker  PR  and  Litvak  MM,  J  Appl  Phys.  1969,  40  {10):  4033. 

[4]  Leiman  HK  and  Oneil  RW,  Appl  Phys  Lett.  1973,  23(1  ):43. 

[5]  David  CS  ,  Proc  IEEE.  1977,  65  {  12  ):  1679. 


24 
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ABSTRACT  Based  on  the  statistic  independent  property  of  the  effects  of  rough  surface  and  at¬ 
mospheric  turbutencc  on  a  propagation  laser  beam,  a  formulation  for  the  probability  density 
lunction  of  the  unwrapped  phase  of  a  elliptic  Gaussian  speckle  field  that  propagates  through 
clear  air  turbulence  is  developed. 

KL\  WORDS  unwrapped  phase,  probability  density  function,  turbulent  atmosphere. 


I.  Introduction 

Multiple  m.odes  [1-3]  have  been  proposed  for  speckle  fields 
formed  by  scattering  of  diffusing  light  beam  particles,  and  phat 
probability  density  function  of  a  light  passing  through  clear 
atm.ospheric  turbulence.  These  m.odels  describe  the  probability 
distribution  of  the  wrapped  phase  fluctuations  of  radian  values 
between  0  and  211  for  the  phase  variation.  However,  in  practical 
applications  of  adaptive  com.pensation  and  precise  ranging,  we 
often  confront  a  situation  in  which  the  phase  fluctuations  of 
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Therefore,  a  study  of 


optical  waves  often  exceeds  2n  radians, 
the  distribution  rule  that  unwrapped  phase  varies  with  surface 
height  fluctuations  and  eddy  current  intensity  when  exceeding  2P. 
radians  for  the  phase  fluctuations  of  light  waves  is  a  rnore 
DTiportant  and  practical  probleia.  This  problem  was  recently 
analysed  by  Draper  [4];  hov7ever,  Draper  et  al.  neglected  the 
noncircularity  of  the  unsm.ooth  and  scattered  echo  waves.  We  knov? 
from,  studying  the  statistical  distribution  theory  and  experim.ents 
subject  to  speckle  fields  of  scattered  echo  waves  from,  unsm.ooth 
boundary  surfaces  that  a  more  objective  distribution  of  the 
speckle  field  of  scattered  echoes  from  unsm.ooth  surfaces  should 
be  elliptical  gaussian  distribution  [5] ,  therefore  the 
ellipticity  param.eter  will  be  introduced  in  the  paper  to  analyze 
the  oscillation  am.plitude  of  the  scattered  echo  wave  field  of  the 
positioning  target,  which  can  satisfy  the  unwrapped  phase 
probability  distribution  after  speckles  (of  elliptic  gaussian 
distribution)  pass  through  the  atm.ospheric  turbulence. 


II.  Theoretical  Analysis 

Since  the  random,  fluctuations  of  the  target  boundary  surface 
and  fluctuations  of  atm.ospheric  turbulence  are  interdependent, 
the  atm.ospheric  turbulence  and  the  target  (unsm.ooth  surface) 
independently  lead  to  phase  fluctuations  of  optical  waves.  In 
analysis,  the  authors  consider  the  interference  between  target 
and  atm.ospheric  turbulence  as  tV7o  independent  disturbance 
factors.  Under  the  situation  of  the  divergent  laser  propagation 
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of  approxiinately— point  sources,  modulation  of  the  transmitted 
optical  beam,  phase  by  the  atm.ospheric  turbulence  can  be 
considered  as  m.odulation  of  the  constant  phase  on  the  target 
scattered  echo  wave.  Thus,  the  unwrapped  phase  probability 
density  function  of  the  scattered  spots  transm.itted  in  the 
atm.ospheric  eddy  current  can  be  derived  from,  the  following 
equation : 


pi  —  )pin)d^i 


(1 ) 


In  the  equation,  Pi4>i p)  is  the  probability  density  function  of  the 
speckle  phase  fluctuation  condition  in  the  absence  of  atmospheric 
turbulence;  p(|j)  is  the  probability  density  function  of  the 
optical  wave  phase  when  the  atm.ospheric  turbulence  are  present. 

Assum.e  that  the  com.pound  vibrational  am.plitude  of  the  echo 
wave  speckle  field  for  target  scattering  is 

£  =  ^ =  /l„ <?"*+/?(/) e +  iv  ( 2  ) 


In  the  equation,  x  and  y  are,  respectively,  the  real  and  the 
virtual  portions  of  the  speckle  field  E;  i|r  is  the  total  phase; 

is  the  constant  vibrational  am.plitude  com.ponent  (not  being 
disturbed)  of  the  light  wave;  and  AV--"  is  the  scattering  random, 
com.ponent  of  the  field.  Since  we  only  can  sim.plify  the  general 
©lliptical  gaussian  distribution  by  coordinate  conversion  into 
the  distribution  product  corresponding  to  interdependent 
com.ponents  [6]  ,  we  consider  only  in  this  paper  the  joint 
probability  density  function  of  the  real  and  the  virtual  portions 
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of  the  following  compound  vibrational  am.plitudes 


pi.\  .,>■)  = 


I  .  Ax“  Av^ 

exp  --^C— 1-  +  -Y-)  =p(x)p{v 

L  O  (J 


Here  ,;^v  =  a-<.v>  is  the  fluctuation  of  x;  <.v>=  A'exp  ( -<r-/2 )  is  the 

inean  value  of  the  x  ccrnponent. 

f;;=<A.v'>=  (..V/2)[l  +  exp( -2fT-’)  -2exp( -c7-)]  is  the 

fluctuation  variance  of  x. 

(T-=  <  Av’>  =  ( A7  2  )  [  1  -  exp  C  -  2(t- ) )  is  the  fluctuation  variance  of  y  [5] 
<...>  indicates  the  statistical  mean  of  the  system.;  N  is  the 
nuinber  of  scattering  eleiTients  of  the  systern;  ^  is  the  ph=^^e 
fluctuation  standard  deviation  of  the  target  scattering  field. 

Eg.  (3)  indicates  the  interdependent  gaussian  distribution 
of  the  X  and  y  com.ponents  of  the  speckles.  This  m.eans  that  Ax 
and  A.v  are  the  dual  random,  circular  com.pound  gaussian  processes 
of  zero  mean  value.  By  introducing  the  ellipticity  parameter 

P  =  2(t;/(ct;+ct2)  ^  relationship  between  the  fluctuation 

intensity  and  of  the  two  com.ponents  and  the  field 
fluctuation  intensity  /=(£-<£>)■  ^  x  .  is 

l  =  T(l  +  p)/2  ,  J  =  7{\-p)I2  ... 


From  probability  theory,  from  the  analysis  in  [6]  and  the 
relationship  7= 2^^^  cos  )  obtained  from.  Eq.  (1)  (as 
shown  xn  Eq.  (5)) ,  we  can  obtain  the  joint  probability  density 
function  between  the  elliptical  gaussian  process  scattering 
intensity  and  the  total  phase  x-^^hen  the  phase  is  constrained  by 
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(5) 


pH  <y)  = 


4n  pa] 


exp 


r  /+/.- 2777*  cosy 

( \  +  p)o] 


—  exp 


i+i-2jTr^  cosy 
{\-p)c] 


In  the  equation  y^>p-n,I==EE'  and /^=  ^tl  =  +  By  using 

the  relationship  between  the  joint  probability  density  and  the 
fringe  probability  density,  we  can  obtain  the  phase  probability 
density  function 


PAy)  = 


3(7;-(t; 

4a:  I  2n 

+  i^i(  ~  7- 


exp 


/,  5in7 

(l  +  p)a;  J 


nl. 


LV  TTTTTTr 


11  Isos-y  ]|  [  1  r  /sin’y 

2  ’  2  ’  (l+p)a;  4cl  [2,1  ""171 -p)a;  _ 


,  cosy-f-,F.(--L--L - )11 

(i-p)a;  ^  2  '  2  ’  (l-p)aj  ^Jj 


(6) 


We  now  consider  the  weak  scattering  target  when  the  surface 
height  fluctuation  is  relatively  small  {  1  );  according  to 

the  surface  light  scattering  model  of  Goodm.an  [7]  ,  the  scattering 
field  phase  i|r  is  proportional  to  the  boundary  fluctuation  height 
tl/  =  a+br}  •  Here  a  and  b  are  constants.  Since  by  now  the  phase 
fluctuation  y~a  +  h)j-fi  is  relatively  sm.all,  then  smy  5:y  ,  cosya  I 
Because  of  a]«\  ,  then  the  confluent  supergeometric  function 

r  1'  _  I  I  ,  ! - 

'I  2  '  ^  V  TT  JC  *  ^  • 

From,  these  relationships,  we  obtain 
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^  (  JL  )  =  1 

'  4a;  L  V  n(l  +  p)a, 


CT*  —  a 


1 


4'^!  L  V  t:  ( 1  —  p  )  a. 


exp(  - 


exp(  - 


y‘ 


(1  +  p 


( 1  -p  )ff 


r  ) 


(7) 


In  the  equation,  o',  =  o]iI,  ■  Actually,  Eq.  (7)  indicates  that 
the  X  and  y  directions  of  the  scattering  echo  waves  of  a  weak 
scattering  target  have  their  phases  of  the  interdependent 
components  satisfying  the  gaussian  distribution.  On  the  other 
hand,  when  the  target  is  a  smooth  plane,  the  surface  fluctuation 
Introduced  by  reflection,  the  additional  phase  is  a 


i=n 


constant.  We  assum.e  that  this  is  zero  and  then  we  have  |j  =  a.  As 
an  approximation,  we  assum.e  a  =ff„  .  Then  from,  the  standard  m.ethod 
[3,4]  to  convert  the  density  function  of  the  unwrapped  phase 
probability  function  into  the  probability  density  function  of  the 
wrapped  phase,  a  sum.m.ation  is  m.ade  on  the  probability  density 
function  pj,  in  the  region  of  (  -  x  .  a)  ) 


/'  (  —  )  =  Z  P*  (  — 


4o' 


\  n  U  *  />  ) 


L  P 


4a; 


V  71  ( I  -  p  )  ap 


I  exp 


{■/+2kn)- 

(1  -  p)ff^ 


(V  2kK  )' 
I  1  +  p  )a,; 


(X  ) 


3y  using  the  sim.ple  harm.onic  decreasing  property  of  exp{~x)‘. 


m.ark  the  equatioi 
for  Eq .  (8) 
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From  studying  the  laser  beam  in  its  transmission  in 
atm.ospheric  turbulence,  when  the  divergent  laser  beam,  is 
transm.itted  and  passes  through  weak  turbulence,  the  instantaneous 
phase  fluctuations  satisfies  the  gaussian  distribution 


p(p) 


(p-Po)~ 

2a‘ 


(10) 


In  the  equation,  <  0^  =  bj  lA'-  ;  b  is  the  random-com.ponent 

variance  of  the  scattering  field  in  the  eddy  current,  and  A'  is 
the  coherent  component  of  constant  vibrational  am.plitude  of  the 
eddy  current  scattering  field.  Under  certain  eddy  current 
conditions,  ol  =  QJii(- LCm^  '  Lj  is  the  exterior  dim.ension  of 

the  turbulence;  L  is  the  distance  from,  the  light  source  to  the 
receiver;  k  =  2nl/.  is  the  num.ber  of  vjaves ;  and  is  the  structural 
constant  of  atm.ospheric  ref ractivity .  Substitute  Eqs .  (9)  and 

(10)  into  Eq.  (1),  and  we  can  obtain  that  the  unwrapped  phase 
probability  density  function  of  the  target  scattering  echo  wave 
in  propagation  in  atm.ospheric  turbulence  is 
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If  the  scattering  echo  wax^e  is 
distribution  speckle,  p=l ,  then  Eq. 
results  obtained  in  reference  [4] . 


circularly  synur.etric  gaussian 
(11)  degrades  into  the 


p('y)  = 


2g' 

(P 


(12) 


Here  G,  =  \a-  +  Gi 


III.  Conclusions 

As  presented  in  the  paper,  the  model  indicates  that  V7hen  the 
phase  fluctuation  of  light  waves  caused  by  turbulence  satisfies 
the  gaussian  distribution  and  the  height  fluctuation  (of  the 
scattering  surface)  is  very  small,  the  phase  distribution  of 
the  unv7rapped  speckle  and  eddy  current  combination,  and  the 
speckle  field  phase  of  the  target  scattering  echo  waves  satisfy 
the  distribution  of  the  same  rule.  This  is  because,  on  the  one 
hand,  the  echo  wave  phase  of  scattering  by  unsm.ooth  boundary 
surfaces  can  be  dissolved  into  the  sum.mation  of  random  phase  and 
constant  phase  satisfying  the  gaussian  distribution.  On  the 
other  hand,  we  can  neglect  the  noncircularity  of  the  speckles  due 
to  the  atmospheric  eddy  current  effect  as  compared  to  the  effect 
of  boundary  surface  scattering.  Therefore,  the  constant  phase  of 


the  boundary  surface  scattering  echo  wave  due  to  the  effect  of 
atmospheric  turbulence  also  becomes  the  variable  of  the  gaussian 
distribution.  Moreover,  since  the  eddy  fluctuations  and  the 
boundary  surface  fluctuations  are  statistically  interdependent, 
V7e  can  obtain  the  above-mentioned  conclusions  from,  probability 
theory.  The  model  also  indicates  when  we  know  the  fluctuation 
variance  of  the  surface  height  of  the  scattering  surface,  the 
phase  variance  of  speckles  or  num.ber  of  scattering  elem.ents  of 
the  target,  and  the  phase  fluctuation  variance  introduce  by  the 
circularity  param.eter  and  turbulence,  we  can  obtain  the  joint 
phase  probability  distribution  rule  of  the  unwrapped  speckle  and 
eddy  current  com.bination .  Eq.  (11)  is  still  adaptable  to  phase 
fluctuation  wrapped  within 

The  first  draft  was  received  on  Novem.ber  3,  1992;  the  final 
revised  draft  was  received  for  publication  on  January  31,  1993. 
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ABSTRACT  The  principle  of  hole  array  wavefront  sensor  is  described  and  the  phase  distribu¬ 
tion  of  a  3,S/im  DF  hii;h  power  laser  beam  has  been  obtained  using  this  method. 

KE^  VVORDS  wavefront  sensor,  DF  laser,  hole  array. 


I.  Introduction 

Because  of  errors  in  machining  the  cavity  lens,  as  well 


c  a  V  i  t  ^ 


thermal  deform.ation  caused  by  laser  energy  absorbed  by  the 
lens,  output  lens,  and  output  window,  and  nonhom.ogeneous 
flowfield,  the  result  will  be  phase  anom.alies  in  the  output  lase] 
beam.  However,  phase  anom.alies  strongly  affect  focusing 
P^^ft^rm. ance  of  laser  beam.s  and  farfield  transm.ission  properties 
[1] .  Phase  inspection  of  lasers  is  one  of  the  key  problems  in 
high-powered  laser  technology.  A  3 . 2-m.icrom,eter  laser  is 
outputted  by  a  DF  chem.ical  laser  device;  the  wavelength  is  at  the 
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mid-infrared  band.  By  citing  an  example  of  the  Hartmann-Shack 
phase  detector  [2],  there  are  m.ainly  three  difficulties: 

(1)  machining  is  required  to  install  the  tuned  infrared  sublens 
array,  at  high  difficulty  and  high  cost.  (2)  A  high  resolving— 
power  infrared  plane  array  optoelectronic  sensor  is  required,  at 
high  cost.  (3)  The  therm.al  effect  of  high  power  will  cause 
additional  phase  anom.alies. 

For  these  reasons,  the  article  adopts  a  sim.ple  hole  array 
phase  sensor  and  presents  its  principle  and  experim.ental  layout 
as  well  as  the  data  processing  method.  The  phase  distribution  of 
the  DF  laser  is  m.easured  directly. 

II.  Principle  and  Experimental  Layout 

The  measurement  principle  is  the  sam.e  as  that  of  Hartm.ann 
wavefront  sensor  [2]  .  The  difference  is  that  the  sm.all  hole 
array  (referred  to  below  as  hole  array)  is  used  to  replace  the 
sublens  array  in  the  Hartm.ann  sensor.  As  shown  in  Fig.  1,  the 
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If  tli0  Isss^ 


incic3,©rit  Isssir  +*1^01  Vi/^io  av'V'OT 

beam  is  a  strict  plane  wave,  then  the  ^y.^^yr  j>i_ 

observation  screen  is  the  same  as  that  of  the  whole  array  (dotted 
lines).  When  the  laser  phase  is  not  planar,  since  the  direction 
of  propagation  of  light  in  very  small  holes  is  perpendicular  to 
the  phase  plane  at  the  site,  the  light  passing  through  +-h‘=“==> 
sm.all  holes  will  propagate  along  solid  lines.  On  the  observation 
screen,  there  will  be  a  translation  of  the  light  spot  array 
position  relative  to  the  planar  wave.  Assum.e  that  the 
translations  corresponding  to  the  i-th  hole  in  the  x  and  y 
directions  are,  respectively,  Ak;,  Ay,.  Moreover,  it  is  assum.ed 
that  the  distance  between  the  hole  array  and  the  observation 
screen  is  f,  then  the  inclination  rate  -  and  g  of  the  wavefront 
of  the  i-th  hole  are: 


g,  = 


A.v 


Ar 


/  f 

therefore,  by  just  .measuring  Ax,  and  Ay,,  and  ,,  can  be 
derived.  Thus,  a  certain  algorithm  can  be  used  to  locate  the 
pha— ^  plane  of  the  laser  at  the  hole  array. 

In  practical  work,  since  there  can  be  two  work  aspects  with 
different  conditions:  (1)  for  infrared  light  with  visible  light 
and  wavelength  less  than  1 . Imicrometers ,  a  silicon  CCD  plane 
array  optoelectronic  sensor  can  be  used  as  the  optoelectronic 
conversion  apparatus  (refer  to  Fig.  2).  The  silicon  CCD  camera 
transform.s  the  light  spot  arrays  on  the  diffuse  reflection  screen 
into  electric  signals,  which  are  converted  to  digital  signals  by 
the  im.age  plate,  then  to  be  inputted  into  a  m.icrocom.puter  for 
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Fig ^2  Case  of  visible  or  infrared  laser  beam 


processing.  (2)  For  a  DF  3 . Smicrorceter  laser  or  a  10 . 6micrometer 
COj,  a  method  similar  to  (1)  can  be  used  when  there  is  a  therm.al 
imaging  instrum.ent  with  corresponding  waveband.  Otherwise,  the 
observation  screen  in  Fig.  1  can  be  changed  into  an  organic  glass 
plate.  When  the  laser  power  is  sufficiently  high,  the  laser  beam 
passing  through  the  hole  array  will  burn  sm.all  holes  on  the 
organic  glass  plate  or  will  leave  visible  molten  tracks.  The 
burn  tracks  on  the  organic  glass  are  inputted  into  the  com.puter 
with  a  scanner.  Method  (2)  is  used  in  this  article. 

III.  Experim.ental  Results 

We  take  into  account  the  spatial  resolving  power  at 
diffraction,  and  the  hole  spacing  of  the  hole  array  is  4 . 0m.m.  and 
the  sm>all  hole  diameter  is  1.5m.m..  To  avoid  the  laser  beam, 
outside  the  hole  being  reflected  into  the  laser  device,  the  sm.all 
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Xins  and  tha  noniTial  Tina  to  tha  hcla  array  plata  ana 
procassad  into  an  angla  of  about  15^  Fig.  3  shows  tha  burn 
track  recorded  on  the  organic  glass. 


Fig.  .1  Burn  hole  on  lucity 


The  spacing  (between  the  organic  glass  plate  and  the  hole 
array)  f=50cm.  The  ratio  between  the  figure  and  the  object  size 
IS  1:1.  The  figure  is  plotted  with  an  HP  Scan jet  Plus  scanner  to 
input  the  information  on  the  organic  glass  plate  into  the 
computer  and  then  vrith  laser  printer  output.  The  resolving  power 
of  the  scanning  instrument  is  300  dots  per  inch;  in  other  words, 
the  spacing  between  the  lines  is  0. 0847m, m.  The  light  spot 
position  on  the  organic  glass  plate  is  given  by  the  center  of  the 
burn  track  shape.  The  wavefront  structure  adopts  the  Southwell 
regional  method  [3].  See  Fig.  4  for  phase  distribution. 


Conclusions 


IV. 


Fig.  4  WaM'I’ront  of  DF  laser  beam  after  reconstruction 

By  using  the  hole  array  method,  the  phase  plane  of  the  DF 
laser  is  directly  measured.  Com.pared  with  the  conventional 
Hartmann  sensor,  the  hole  array  method  has  the  following 
advantages:  (1)  the  equipm.ent  is  quite  sim.ple.  It  only  requires 
a  Kietal  plate  and  an  organic  glass  plate  to  be  drilled  with  a 
small  hole  array,  without  requiring  precise  adjustment.  The  cost 
is  very  lov7.  (2)  The  hole  array  can  be  used  in  visible  light  or 
in  infrared  light.  To  a  certain  extent  the  difficulties  of  phase 
m.easurem.ent  of  infrared  lasers  are  solved.  (3)  High-powered 
lasers  can  be  directly  m.easured  without  deterioration.  (4)  There 
is  higher  spatial  resolving  power.  (5)  In  the  optical  path,  no 
additional  optical  elem.ents  are  required,  therefore  there  is  no 
additional  phase  anom.aly.  However,  there  are  som.e  problem.s 
remaining  to  be  im.proved:  one  is  the  precise  measurem.ent  of  the 
light  spot  position.  Secondly,  the  tim.e  dependent  variation 
difficulties  are  greater  for  m.easuring  the  wavefront  by  using  the 
burn  track  m.ethod  with  an  organic  glass  plate.  Equipm.ent  sim*ilar 
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to  that  of  Fig.  2  can  replace  the  organic  glass  burn  hole  V7ith  a 
therrr.al  instruinent.  The  authors  are  developing  studies  in  this 
area . 

The  first  draft  of  the  paper  was  received  on  January  31, 
1993,  The  final  revised  draft  was  recei^’’^'^  ^"'^'^il  1  iqq7 
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ABSTRACT  The  recently  appeared  several  neve!  optical  unstable  resonators  are  reviewed.  Their 
characteristics  and  development  arc  analysed. 

KEY  WORDS  novel  optical  unstable  resonators. 

I.  Introduction 

For  large-energy  high-powered  laser  devices,  the  most 
important  matter  is  dealing  with  the  contradiction  between  laser 
output  energy  and  light  beam  quality.  The  optical  resonator  in 
practice  follows  tradeoff  in  selective  principles  between  energy 
output  and  light  beam  quality  to  achieve  large-mode  with 
prerequisite  light  beam  quality  of  indicators  to  obtain  near¬ 
diffraction  limit  Gaussian  beam  operating  at  TEMqq  mode  in  the 
small  nearby  region  of  the  optical-axis  center  for  the  classical 
stabilized  cavity  in  intermediate  and  low-powered  laser  devices. 
However,  for  large  mode  size,  in  addition  to  prolonging  the 
cavity  length  but  also  using  an  in-cavity  optical  beam-widening 
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element,  this  will  induce  phase  anomalies  due  to  image  difference 
of  beam  expansion  elements  and  optical  element  easily  damaged  by 
the  laser  in  the  cavity. 

However,  for  laser  devices  with  large  Fennell  number  it  is 
not  practical  to  obtain  a  large  volume  output  of  the  fundamental 
mode.  For  the  ring-shaped  outputted  concave  and  convex 
stabilization  cavity,  if  the  parameter  g  can  be  strictly 
controlled  within  the  criterion  range  of  the  stabilization 
region,  theoretically  speaking  a  greater-size  fundamental-mode 
output  can  also  be  obtained.  Since  the  design  conditions  are  too 
strict,  only  when  A=L+Rcoiivex"’^concave  approaches  zero,  and  gigj— >1, 
can  the  above-mentioned  output  be  obtained  at  the  fringe  of  the 
stabilization  region.  Therefore,  the  machining  precision  and 
alignment  adjustment  of  the  resonator  are  especially  difficult. 
Later,  such  cavities  were  replaced  with  conventional  confocal 
unstable  cavities , 

The  conventional  confocal  unstable  cavity  can  be  used  to 
obtain  large-mode  volume;  however,  the  ring-shaped  output  beam  is 
quite  sensitive  to  maladjustment  of  the  mirror  and  anomalies  in 
the  mirror  surface,  thus  causing  inhomogeneity  of  the  laser 
patterns  in  the  near-field  of  the  ring-shaped  beam.  It  is  very 
damage  the  window  of  a  high-powered  laser  device. 

However,  the  far-field  light  beam  patterns  still  unavoidably 
produce  many  diffraction  rings  affecting  light  beam  quality.  On 
account  of  this,  the  divergence  of  the  10,000-W  class  CW  CO, 
laser  device  by  using  the  conventional  unstable  cavity  is  still 
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unable  to  be  attained  within  several  times  the  diffraction  limit. 
This  ring-shaped  diffraction  output  cavity  has  another  flaw. 
Within  the  range  of  the  diffraction  bandwidth  (A./©)  at  the  fringe 
of  the  mirror  (©  is  the  diffraction  output  angle  of  the  light 
beam  at  the  fringe  of  the  mirror) ,  since  the  thermal  deformation 
of  the  lens  or  due  to  processing  and  installation,  the  angle  ©  is 
sufficiently  large,  some  light  beams  outputted  with  oscillation 
diffraction  of  the  laser  in  the  cavity  will  return  to  the  center 
zone  of  the  optical  axis  in  the  cavity  along  the  route  of  its 
coming,  to  excite  the  new  low-valency  transverse  mode  period, 
inducing  a  complicated  mode  structure.  This  situation  will  be 
more  serious  with  an  increase  in  laser  power,  in  particular,  for 
lower-gain  laser  devices  adopting  the  low  amplifying-rate  design 
to  upgrade  the  light  beam  quality  and  eliminate  the  side  lobe 
effect,  as  well  as  to  eliminate  the  effect  on  complex  mode 
structure  in  the  cavity  due  to  the  hard  edges  of  the  lens.  Thus, 
the  various  low-order  modes  separate  and  achieve  the  complete 
fundamental  mode  operation  require  the  development  of  novel 
unstable  cavity  with  large-mode  volume. 

In  recent  years,  on  several  occasions  abroad  proposals  for 
unstable  cavity  methods  of  large  mode  volume  to  greatly  upgrade 
the  light  beam  quality  of  high-powered  lasers  were  proposed. 

This  article  will  detail  and  discuss  these  novel  unstable 
cavities  with  application  prospects  given. 
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II.  Several  Novel  Unstable  Cavities 

2.1.  Radial  variable— ref lectivity  mode  unstable  resonator 
(RVRMUR) 

This  cavity  is  also  called  the  unstable  cavity  with  a 
Gaussian  reflectivity  output  coupler  (GRC) .  This  cavity  adopts 
the  step-by-step  output  approach  of  transmission.  Therefore,  the 
effect  on  the  fundamental  mode  operation  by  complex  mode 
structure  in  the  cavity  due  to  the  hard  edges  of  the  conventional 
unstable  cavity  can  be  eliminated.  In  the  method,  by  adjusting 
the  reflectivity  at  the  mirror  surface,  the  distribution  of 
mirror  surface  reflectivity  is  of  a  Gaussian  function.  In  1970, 
H.  Zucker  et  al.  were  the  first  to  propose  this  concept  [1]. 
However,  since  the  output  mirror  technique  has  not  been 
qualified,  this  technique  was  not  emphasized  in  the  academic 
field.  Only  several  researchers  conducting  experiments  confirm 
this  theory  on  laser  devices  with  relatively  low  power  and 
relatively  small  cross-sectional  area  [2-5] . 

In  1985,  P.  Lacigne  et  al .  in  Canada  coated  the  high 
reflectivity  medium  mode  varying  with  the  mirror  diameter  onto 
the  substrate  with  an  antireflection  treatment.  In  addition, 
denaturation  with  transmission  of  nonreflected  light  to  remedy 
the  unstable  cavity  that  the  Gaussian  reflectivity  mirror  is  used 
as  the  output  coupler.  Then  this  technique  was  applied  in  high- 
power  CO2  laser  devices.  By  applying  the  output  mirror  with  the 
Gaussian  reflecting  surface,  the  harmful  effect  of  mirror  edge 
diffraction  in  conventional  unstable  cavity  can  be  reduced.  In 
addition  the  second  diffraction  ring  in  the  distribution  of  far- 
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field  intensity  is  eliminated,  thus,  upgrading  the  directivity 
and  brightness  of  the  light  beam  [7].  In  1988,  S.  deSilvestri  et 
al .  proposed  a  super-Gaussian-ref lectivity  mirror  [8].  The 
fundamental-mode  volume  of  the  unstable  cavity  with  this  output 
mirror  is  greater  than  the  former  one,  and  the  laser  output 
efficiency  is  also  higher.  In  a  high-Fennell-number  resonator, 
good  light  beam  quality  can  also  be  obtained  although  from  the 
super-Gaussian  reflective  type  surface  equation,  we  obtain 

/?,=,R,exp[-2(^  )’]  (n>2) 

m  ^ 

In  the  equation,  Rq  is  the  reflectivity  at  the  center  of  the 
output  coupler;  r  is  the  radius  of  the  mirror;  and  Wjj  is  the 
radius  at  the  GRC  reflectivity  curve  1/e.  We  can  see  that,  with 
increase  in  n,  the  super-Gaussian  indicator,  the  diametral 
direction  reflectivity  type  surface  curves  gradually  have  the 
flat-topped  phenomenon,  which  is  disadvantageous  to  mode 
selection.  However,  even  if  9  is  selected  for  m,  in  the  COj  and 
YAG  laser  devices,  experiments  proved  that  a  super-Gaussian 
mirror  can  still  produce  beams  with  near-diffraction  limit,  and 
can  output  a  higher  energy  over  the  classic  unstable  cavity. 

This  cavity  is  composed  of  a  piece  of  concave  surface 
100%— ref lecti ve  mirror  and  a  piece  of  output  coupling  mirror  with 
variable  reflectivity  output  at  the  diametral  direction  of  the 
convex  surface.  In  this  cavity  type,  light  beam  waist  in  the 
cavity  is  determined  by  the  following  equation: 

=  (2) 
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In  the  equation,  M  is  the  amplifying  rate  of  the  optical  cavity; 
Wjj  is  the  radius  of  the  GRC  reflectivity  curve  at  1/e. 


Fig.  I  Uns(ah!o  resonator  wuh  eaussian 
rcdcctivity  mirror 

By  varying  the  mirror  waist  W,,  and  the  optical  cavity 
amplifying  rate  M,  and  by  using  Eq.  (1),  we  can  obtain  any 
Gaussian  beams.  The  distribution  of  the  near-field  optical-beam 
intensity  after  using  this  cavity  type  can  be  expressed  by  using 
the  difference  of  two  Gaussian  functions: 

A.o=--/,,<,exp{  :[1 -R„exp'l  -  }]  (3) 

In  the  equation,  r  is  the  diametral  distance  from  the  optical 

axis,  and  is  the  light  intensity  at  the  GRC  center  optical 

axis  . 

If  the  center  reflectivity  /?^>  ( 1  +  iC-  )-'  sinking  at  the 
center  will  appear  in  the  near-field  distribution  (Fig.  2).  The 
dotted  lines  in  the  figure  indicate  the  unstable  cavity,  while 
solid  lines  indicate  the  unstable  cavity  with  a  phase-unifying 
coupler.  When  the  transmission  beam  is  straightened  with  the  GRC 
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and  the  far-field  brightness,  when  1  ,  is  the  highest,  when 

far-field  intensity  can  be  approximately 
determined  by  the  following  formula 


r.9  ^/.o  ' 


exp{  - 


er 


2(1  + 


jyr-«p 


_ A 

Wi 


{■  L 

^2  '  ^  0 


(4) 


In  the  equation,  &a(—^-l2nW^)  is  the  divergent  half-angle  at  the 
maximum  intensity  1/e  of  the  Gaussian  beam  in  the  cavity;  If ^  is 
the  light  intensity  at  the  center  optical  axis  of  the  far  field. 

From  the  far-field  intensity  distribution  diagram,  we  can 
see  that  the  second  diffraction  ring  of  the  far-field  light  beam 
has  been  severed  (Fig.  3).  Another  advantage  of  the  unstable 
cavity  of  this  diametral  variable-reflectivity  mirror  is  its 
nonrestraint  by  the  gain  cross-section.  Although  the  mismatching 
of  the  cavity  adjustment  will  lower  the  output  energy  of  the 
light  beam,  but  it  does  not  affect  the  output  of  the  far-field 
solid-core  Gaussian  type  surface. 

Fig.  4  shows  the  fundamental  design  of  the  GRC  reflected- 
output  mirror.  The  output  mirror  is  composed  mainly  of 
multilayer  high-ref ractivity  medium  film  with  a  certain  shape 
chemically  deposited  on  the  anti-reflection  medium  substrate. 

Some  of  the  unreflected  light  passes  through  the  anti-reflecting 
substrate.  Along  with  the  diametral  variation,  the  thickness  of 
the  medium  film  surface  can  be  derived  from  the  theory  of 
multiple  light-beam  interference  of  the  thin  film  of  the  medium. 
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Fig.  2  (A)  Intracavity  and  ( B )  near- field  intensity  distributions  of  the  beam  produced  in 
a  resonator  ( A/=  2 )  with  a  GRC  ( 1 .2cm»  =  0.7 ) 


4nn2 


\-R. 


-’"nr.,- 


(5) 


In  the  equation  Rj,  is  the  diametral  reflectivity;  is  the 
designed  wavelength  in  vacuo;  n2  is  the  reflectivity  of  the 
medium  film;  r|2  is  the  reflectivity  of  the  surface  amplitude  in 
the  medium  between  1  and  2;  rg^  is  the  amplitude  reflectivity  at 
the  inner  surface  of  the  antireflecting  film  substrate. 

Generally  in  abroad,  the  film  coating  of  such  mirrors 
applies  the  chemical  vapor  deposition  (CVD)  method.  In  other 
words,  this  is  the  application  of  a  gaseous  substance  to  have  a 
chemical  reaction  at  the  solid  surface  in  the  process  of  forming 
a  solid  deposit.  Therefore,  the  affinity  between  the  formed  film 
and  the  substrate  is  better  than  the  physical  vapor  deposition 
method  (PVD) .  The  CVD  process  can  take  place  at  constant 
pressure,  and  also  at  low  pressure.  Since  at  low  pressure  the 
diffusion  rate  of  a  gaseous  reflecting  substance  is  greater  than 
that  at  constant  pressure,  this  has  a  good  function  on  the 
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Fig.  3  Far -field  intensity  distribution  of  the  beam  coupled  out  of  a  resonator  (A/  = 
i:.5cm- diameter  GRC  (2 ff'.=  2.4cm ) and  (B)  around  a  2.5cm -diameter  hard 


2 ) ;  ( A )  through 
“edged  mirror 


a 


(n,) 


PreflUd  lajrvr  of  ditleetrlc 


>n-cooiin„ 
inj) 

Fig.  4  Configuration  of  the  complementary  Gaussian  rcficctivity  mir 


diffusion,  adherence,  reaction,  and  growth  of  film  substance 
between  the  reactive  substance  and  the  matrix.  Therefore,  there 
is  obvious  improvement  in  the  quality  of  the  film  layer.  Thus, 
in  recent  years,  the  low-pressure  method  was  applied  in  CVD. 
Generally,  using  CVD  in  making  thin  films  is  divided  into  four 
main  stages:  (1)  surface  diffusion  of  the  reactive  gas  toward  the 
matrix;  (2)  adherence  of  the  reaction  gas  onto  the  surface  of  the 
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matrix;  (3)  chem.ical  reactions  occur  at  the  matrix  surface; 

(4)  on  the  surface  of  the  matrix,  the  gaseous  by-products 
produced  are  detached  from  the  surface  and  are  pumped  out  and 
discharged  by  the  vacuum  pump,  so  a  firm  reactive  substance 
remains  on  the  substrate  surface.  Previously,  the  CVD  had  not 
been  extensively  promoted.  The  main  reason  was  the  higher 
deposition  temperature.  In  recent  years,  many  specialists 
conducted  studies  on  lowering  the  thin  film  deposition 
temperature  and  explored  new  reaction  excitation  methods.  Good 
progress  was  obtained.  The  more  successful  of  these  methods  is 
the  molecular  beam  extension  technique,  combining  the  plasma 
excitation  chemical  vapor  deposition  method  (FECVD) ,  and  the 
laser  induced-chemical  vapor  deposition  method.  Therefore,  in 
present-day  international  practice  the  CVD  method  has  replaced 
the  PVD  method  in  making  high-energy  and  high-powered  cavity 
mirrors.  In  the  process  of  CVD  preparing  of  the  diametral 
variable-reflectivity  cavity  mirror,  the  key  technique  is  that  of 
(1)  for  the  development  and  determination  of  the  dynamic 
parameters  of  the  chemical  gaseous  source  in  the  components  of 
the  coated-film  system  and  the  reaction  process  (the 
concentration,  reaction  temperature,  and  pressure  of  the  chemical 
reactant).  (2)  Study  and  determine  the  shape  parameters  of  the 
shielding  apparatus  in  the  CVD  chamber,  and  the  relative  position 
of  the  deposited  substrate. 

From  the  foregoing,  by  using  the  diametral  variable- 
reflectivity  unstable  cavity  coupled  with  the  GRC,  the  side  lobe 
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pattern  in  the  far-field  intensity  distribution  can  be  eliminated 
in  order  to  upgrade  the  light  beam  directivity  and  brightness. 
This  is  attractive  for  the  low-gain  and  low  amplif ying-rate 
resonators  with  the  second  diffraction  ring  occupying  the  higher 
proportion.  The  machining  of  this  type  of  cavity  mirror  is 
complex.  Currently,  there  are  only  a  few  countries  in  the  world 
that  have  developed  successfully.  In  particular,  there  have  been 
no  development  reports  on  the  variable-reflectivity  mirror  with 
X-Y  axial  symmetry  with  rectangular  output  light  beam  and  light 
cavity.  In  the  author's  view,  although  this  film  coating  machine 
is  more  complicated  in  structure  and  processing  technique,  yet 
this  is  a  new  approach.  Its  success  will  mark  the  occasion  of  a 
high  degree  of  unification  of  light  beam  quality  and  laser  output 
efficiency  of  the  intermediate  and  low-gain  transverse-flow  laser 
devices . 

2.2.  Phase-Unifying  Output  Coupler  with  Unstable  Resonator 
(PUCUR) 

In  1988,  S.  Yagi  et  al.  in  Japan  were  the  first  to  publish  a 
paper  on  an  unstable  resonator  with  phase-unifying  output 
coupler .  This  is  also  called  an  unstable  resonator  with 
diametral  graduated  variable-reflectivity  mirror  [9] .  Such 
cavity  has  partial  reflection  from  its  center,  composed  of 
output-coupled  mirror  and  an  all-concave  reflecting  mirror  with 
anti-reflectivity  at  the  fringe  (Fig.  5) .  The  film-coating 
processing  method  for  this  output  mirror  is  actually  a  variety 
and  a  simplification  of  the  diametral  variable  super-Gaussian- 
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reflectivity  output  mirror. 


Its  selective  capability  in  the 


Fig-  5  Unstable  resonator  with  phase- unifying  coupler 


transverse  mode  is  less  than  that  of  the  unstable  resonator  of 
the  super-Gaussian-ref lectivity  mirror  with  smaller  m.  However, 
its  mode  volume  utilization  is  close  to  that  of  the  resonator 
with  a  greater  super-Gaussian  reflectivity  mirror.  So  such 

^®sonator  is  more  suitable  to  the  CO2  laser  device  with 
longer  wavelength  and  intermediate  gain.  Since  this  also  adopts 
the  output  approach  with  transmittance  or  m.ainly  with 
transm.ittance ,  therefore  this  does  not  have  the  harmful  effect  of 
diffraction  of  the  conventional  mirror  with  unstable 
resonator.  Such  unstable  resonator  can  be  used  to  induce  the 
phase-unified  far-field  light  beam  with  large  volume  of  desirable 
pattern . 

On  the  initial  field  distribution  of  the  diametral  light 
beam,  on  the  surface  of  this  type  of  resonator  mirror,  let  us 
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assume  that  the  pseudorandom  distribution  of  the  field  vibration 
amplitude  and  phase  on  this  mirror  surface  is  U(rj)  ,  with  Fennell 


(a)  (b)  (c) 

Fig.  6  Comparison  of  laser  power  and  far -field  beam  pattern 
(a)  laser  power  vs  discharge  power,  and  far -field  patten  of  (b)ihe  novel  resonator,  and  (c) 
Conventional  unstable  resonator  at  a  discharge  power  of  43k\V 


integration,  the  field  distribution  U(r2)  on  the  mirror  of  the 
opposite  surface  can  be  determined*  After  the  light  beam  makes 
several  round  trips,  the  field  distribution  in  the  resonator 
reaches  the  saturation  state,  then 


t/(rj=27r 


\r,\<a,Kir,,r,)U(r,)r,dr, 


but 


K 


( r, ,  r, )  =  VU,  (r,  ) (r,)]L 


(6) 

(7) 
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In  the  equation  rj  and  r2  are  the  diametral  distances  of  the 
mirror;  L  is  the  resonator  length;  k=X/2n,  X  is  the  wavelength; 
aj  is  the  radius  of  the  reflected  disk;  and  functions  g^ (r)  and 
g2(r)  are  the  gains  of  the  oscillation  amplitude  in  front  of  each 
mirror . 

Assume  that  the  saturation  gain  is  given  by  the  formula: 

_ |o _ 

In  the  equation,  gj)  is  the  small-signal  gain  of  the  medium;  I 

3 

and  are  the  local  field  intensities  of  the  plane  wave 
propagation;  is  the  saturation  strength. 

The  complex  oscillation  amplitude  function  relates  to  the 
electric  field 

£(r)  ~Rct/(r)cxp( 

In  the  equation,  Wp=2nc/Aa,  c  is  the  speed  of  light. 

In  the  resonator,  intensity  I(r)  and  phase  P(r)  can  be 
calculated  from  the  following  equations: 


/(r)  =  |{/(r)P 

/’(/•)  =  arctan{  — } 
RcU(r)  ‘ 


( 10) 
(II) 


Fig.  5  shows  the  light  beam  patterns  of  the  near-  and  far- 
fields  after  the  CO2  laser  adopts  such  kind  of  resonator.  By 
m.aking  a  comparative  experiment,  it  was  discovered  that  by  using 
this  new  unstable  resonator,  its  power  output  is  slightly  less 
than  that  of  a  conventional  unstable  resonator.  In  the  former, 
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the  power  output  level  of  M=1 . 5  approaches  that  of  a  conventional 
unstable  resonator  designed  on  the  basis  of  M=2.  However,  the 
high  focusing  beam  of  the  solid  pattern  generated  by  this  cavity 
(M=1.5)  is  better  than  that  of  a  conventional  resonator  {M=3) 

(fig.  6)  [10].  S.  Yagi  et  al .  of  the  Mitsubishi  Electric 

Corporation  in  Japan  obtained  linearly-polarized  light  with  a 
divergence  angle  of  0 . 6mrad  with  a  5kW  CO2  coronal  discharge 
transverse-flow  laser  device.  Based  on  the  output  aperture  of 
45mm,  calculation  yielded  the  diffraction  limit. 

2.3.  Self-filtering  Unstable  Resonator  (SFUR) 

By  adopting  the  high-amplif ying-rate  design  for  conventional 
high-gain,  high-power  laser  device,  its  conventional  unstable 
resonator  can  produce  better  light  beam  quality.  However,  this 
is  very  difficult  to  adopt  the  high-quality  laser  beam  by  using 
the  same  approach  for  the  high-power  short-pulse  excimer  laser 
device.  The  reason  is  that  the  gain  buildup  time  for  such  laser 
device  is  considerably  shorter  than  the  mode  buildup  time, 
therefore  it  is  unable  to  reach  the  diffraction  limit. 

In  1984,  P.  G.  Cobbi  et  al.  in  Italy  [11]  were  the  first  to 
realize  a  negative-branch  unstable  resonator  of  the  self¬ 
filtering  unstable  type  (SFUR)  (Fig.  7).  Such  resonator  includes 
a  field-limiting  aperture  (LA) .  The  small-hole  radius  (of  the 
aperture)  q  5 jy- /  •  The  resonator  is  situated  at  the 

confocal  point  with  focal  lengths  fj  and  fj  of  the  confocal 
points  of  two  concave-surface  reflecting  mirror  Mj  and  Mj .  When 
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the  plane  wave  impinges  on  LA  and  reflects  back  from  the  mirror 


Fig.  7  Schemes  of  unstable  resonators 

(a  )self- imaging  unstable  resonator :  (b)modilled  form  of  self- filtering  unstable  resonator 


with  smaller  radius  of  curvature,  only  the  Airy  disk  is 
allowed  to  be  propagated  outside  of  the  hole,  and  makes  a  round- 
trip  inside  the  resonator  with  amplifying  rate  M=-f2/f2.  This 
scheme  is  especially  suitable  to  high— gain  laser  devices  such 
that  the  stable  low-order  mode  can  be  quickly  established  in  the 
resonator.  However ,  when  SPUR  is  used  for  a  short— pulse  excimer 
laser  device,  generally  the  mode  volume  is  very  small.  To  let  it 
fill  the  entire  active  volume  cross-section  in  order  to  obtain 
high-level  energy,  recently  P.  D.  Lazzaro  [12]  proposed  a 
generalized  self-f iltering  unstable  resonator  (Fig.  8). 

To  obtain  the  SPUR  property  with  the  same  feature,  it  is 
required  to  introduce  two  conditions:  (1)  self-imaging  condition; 
this  is  the  relationship  between  1|  and  Ij  ,  on  the  one  hand,  and 
f2,  on  the  other,  in  order  to  satisfy  the  requirement  that  the 
mirror  image  of  the  needle  hole  corresponds  to  the  surface  of  the 
needle  hole  itself  after  one  round  trip  between  the  mirror 
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surfaces.  (2)  To  interrupt  the  oscillation  of  the  TEM, 


Fig.  8  Schematic  diagram  of  the  GSFUR 
.U,  ;  aluminized  flat  mirror;  Af, :  aluminized  concave  mirror  of  focal  length /j; 

OC  :  output  coupler ;  LA  :  field  -  limiting  aperture  ;  Af  j :  aluminized  collimating  mirror  of 
focal  length/,;  distance  of  A/,  from  OC 


intrinsic  mode  of  the  shorter  portion  of  the  light  path  separated 
by  the  LA  in  the  optical  cavity,  it  is  required  to  select  the 
needle  hole  radius. 

If  L,f,(2g,-1)= (12) 


The  condition  of  (1)  can  satisfy  the  requirement  that  Lj  and 
in  the  equation  are,  respectively,  the  distances  of  and  M2 
from  the  aperture  LA. 


,  -i-A. 

*’  2/, 


■  *  V, 


(13) 


For  the  planar  mirror  Mj,  gi=l.  From  Eq.  (12),  we  can  obtain 


(14) 


To  interrupt  the  oscillation  of  the  intrinsic  mode  oscillation 
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(TEM(,(., )  in  the  Lj  region,  the  radius  of  the  LA  aperture  should  be 

°  =  (15) 

In  this  cavity,  the  amplifying  rate  is 

. - ^ 

(16) 

From  the  above  equation,  we  can  see  that  the  amplifying  rate 
can  be  varied  by  changing  the  distance  to  Lj .  From  Fig.  8,  we 
can  see  that  there  is  some  bending  of  the  wave  array  surface  of 
modes  in  the  light  cavity.  To  obtain  the  output  light  beam  near 
the  needle-hole  screen,  the  alignment  output  should  require  the 
reflective  mirror  ,  with  focal  length  as 


2g,L, 

2g.-\ 


(17) 


In  the  equation,  d  is  the  distance  between  aperture  LA  and  M-, . 

For  this  self-filtering  unstable  resonator  in  the  broad 
sense,  to  obtain  large  mode  volum.e,  the  am.plifying  rats  should  be 
very  large  (generally,  M  is  greater  than  or  equal  to  10) .  In 
such  unstable  resonators  have  been  successfully 
operating  in  the  Nd:YAG  lasers  and  the  TEA  COj  excimsr  lasers. 

In  particular,  the  unique  feature  can  be  better  exploited  in  the 
high-gain,  high-power,  short-pulse  excim.er  device  that  the  stable 
low-order  m,ode  is  unable  to  be  established  quickly.  After  the 
ultraviolet  pre-ionized  XeCl  laser  device  adopts  the  plane  and 
concave  self-filtering  unstable  resonator,  the  diffraction-lim.it 
Gaussian  laser  beam,  has  been  obtained  for  11ns  short-pulse 


lasers,  V7ith  brightness  at  4.8x10*"  Wcin'^sr**.  If  wave  chopping  or 
other  technique  of  complex  pulse  width  is  adopted  to  increase  the 
brightness  of  the  laser  apparatus  by  two  orders  of  magnitude, 
this  is  very  possible  to  be  used  in  research  on  ICF.  However,  it 
is  required  to  have  high  m.agnifying  power  for  this  type  of 
cavity,  and  this  is  not  suitable  for  low-gain,  high-power  lasers. 
This  is  also  of  the  negative-branch  unstable  resonator  and 
needle-hole  surface  plate  in  the  cavity  that  is  very  easily 
dam.aged.  Due  to  the  cavity  length  being  long  in  the  experim.ents , 
requirem.ents  on  straightening  of  the  cavity  are  also  quite  high. 


Table  I  Advances  in  unstable  resonators 


No 

resonator  type 

laser  set 

outpower 

Or: 

1  1 

GRC  VRUR 

CW  CO, 

300W 

0.35mrad 

2 

GRC  VRUR 

pulse  YAG 

385mJ,  0.5  ^5Hz 

0.26  —0.69mra<l 

3 

PUCUR 

CW  CO, 

5kW 

0.6mrad  <PA5 

4 

GFSUR 

i 

short”pulse 

XcCI 

3.5mJ  1 1  ns 

4.8  X  Ky’Wcm-'sr*' 

0,26mrad 

III.  Conclusions 

The  three  above-m.entioned  new  models  of  unstable  resonators 
have  good  prospects  in  the  high-energy  (or  high-power)  large- 
m.ode— volume  unstable  resonators .  Table  1  shows  the  developm.ent 
level  obtained  at  present. 

In  Table  1,  the  GRC  VRUR  is  suitable  to  high-powered  lasers 


59 


of  intermediate  and  low  gain.  The  advantage  can  be  more 
exhibited  in  the  low^-gain  laser  device.  PUCUR  is  suitable  to  the 
interm.ediate-gain  laser  devices.  SPUR  and  GSFUR  are  suitable  for 
high-gain  laser  devices.  Other  than  the  three  above-mentioned 
new  unstable  resonators,  recently  there  has  been  the  beam- 
rotating  laser  device  and  the  diffraction  filtering  unstable 
resonators,  represented  by  UR9Q.  It  is  noteworthy  that  vigorous 
research  is  also  underway  in  the  United  States. 

The  first  draft  was  received  on  Decem.ber  21,  1992;  the  final 
re\'ised  draft  was  received  for  publication  on  April  5,  1993. 
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